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Abstract 

Properties of transport of molecular motors are investigated. A simplified model based on the concept of Brownian ratchets is 
D applied. We analyze a stochastic equation of motion by means of numerical methods. The transport is systematically studied with 
£h respect to its energetic efficiency and quality expressed by an effective diffusion coefficient. We demonstrate the role of friction 
I and non-equilibrium driving on the transport quantifiers and identify regions of a parameter space where motors are optimally 
i-^ transported. 
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1. Introduction 

In the process of evolution Nature has created perfect 
biomachines. Researchers (read: visionaries) have dreamed 
to get ahead of Nature and construct molecular-level ma- 
chines which would be at the ultimate limit of miniatur- 
ization. Recent development of nanofabrication technol- 
ogy establishes a foundation for designing, synthesizing, 
constructing and testing functional hybrid mechanical and 
electrical devices on a nanometre scale. Such nanomachines, 
operating under different conditions, could find very many 
new applications. Nowadays, there are many problems that 
have to be overcome: how to power such systems, actuate 
and transport them, and couple to another systems. Na- 
ture prompts and provides partial solutions: Within every 
living cell there is a service structure for its effective func- 
tion and operation, transport and distribution of various 
products and substances inside, outside and across cells 
(Howard, 2001). These processes are mediated by tiny bio- 
motors that move along filamentous tracks. They do me- 
chanical work consuming energy taken from the hydroliza- 
tion of adenosine triphosphate (ATP) which operates as a 
power station. The supcrfamilies of molecular bio-motors 
like kinesin, dynein and myosin are responsible for trans- 
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port of vesicles and organelles, chromosome segregation, 
cell division and muscle contraction, to mention only a few. 
Although biological, chemical and mechanical aspects have 
been studied, relatively little is known about the mecha- 
nisms how the biological motors work. The fundamental 
question is whether their movement is purely deterministic 
or rather random which only looks like deterministic on a 
macro-scale. 

Biological motors possess many of the characteristics 
required to power molecular machines. They can generate 
force and torque, transport various cargoes and are able 
to operate in a processive manner, i. e. they can move 
continuously along the specific substrates for distances of 
up to hundreds of steps (several microns). An example of 
such a motor is conventional kinesin (Howard, 2001). The 
two heads of the kinesin dimer work in a coordinated man- 
ner to move along one of 13 protofilament tracks of the 
microtubule. Each protofilament consists of asymmetric 
a/3-tubulin heterodimers. A heterodimer is about 8 nm 
long and is asymmetric because it is composed of two glob- 
ular subunits: a-tubulin and /3-tubulin which are joined 
together in a head-to-tail fashion so that the dimmer has 
a translational symmetry. Because the a/3-tubulin het- 
erodimers are asymmetrical, the microtubule is polar and 
its ends are structurally different. One consequence of this 
polarity is that polymerization is faster at one end than 
the other. The fast-growing end is called the plus-end, 
whereas the slow-growing end is called the minus-end. The 
conventional kinesin moves towards the plus-end. There 
are other bio-motors like e.g. ncd dimer (Henningsen and 



Schliwa, 1997; Sablin et al., 1998) which move in the oppo- 
site direction on the same structure. The reason, why the 
motor moves in a given direction, is not explained to the 
end. There are several conjectures supported by analysis 
of very simple mechanical models. One of them says that 
it can depend on size of the motor, another says that it 
can depend on mass of the motor or more enigmatically it 
depends on " molecular topology of the motor domain rela- 
tive to the rest of the molecule" (Hcnningsen and Schliwa, 
1997). 

Artificial molecular motors are under constructions and 
experiments on transport properties are performed. An ex- 
ample is a Brownian motor by using cold atoms in a dissi- 
pative optical lattice. In experiments the optical potential 
is spatially symmetric and the time symmetry of the sys- 
tem is broken by applying appropriate zero-mean ac forces. 
As a result, a current of atoms through the optical lattice 
can be generated (Jones at al., 2004; Sjolund et al., 2006). 

2. Model 

From a physical point of view, the kinesin (or other bio- 
motors) moves in a one dimensional spatially periodic po- 
tential V(x) = V(x + L) of period L rts 8nm (Visscher 
et al., 1999) . The reflection symmetry of this potential is 
broken V(x) ^ V(—x), because the ab-heterodimers, that 
build the microtubule, a highway for the biomachine, are 
asymmetrical. This is a crucial fact because it determines 
the mechanism of the kinesin movement the ratchet effect 
(Astumian and Hanggi, 2002; Hanggi et al., 2005; Luczka 
et al., 1995). To explain this mechanism, let us consider as 
simple model as possible. We assume that the molecular 
motor is a particle of mass m moving in a periodic potential 
V(x) of period L and of the barrier height AV = V max — 
V m i n . The equation of motion for the motor is the Newton 
equation with a complementary random force which corre- 
sponds to thermal fluctuations, i.e. the Langcvin equation 
in the form 



mi + 7± = f(x) + g(t) + V2DT(t) 



(1) 



Inertial effects, related to the mass, are described by the 
first term on the left hand side. The dissipation is included 
via the Stokes force with the friction coefficient 7 which is 
proportional to linear size R of the particle i.e., 



7 = 6irr)R 



(2) 



and is additionally determined by the viscosity r\ of the 
medium the particle moves in. The potential force 



m = --^ (3) 



dV(x) 
dx 

is zero over a period L, 

L 

(f(x)) L = ± f f(x)dx=±[V(L)-V(0)}=0. (4) 



The stochastic force T(t) describes thermal fluctuations 
which can be modeled by (5-correlated Gaussian white noise 
of the statistics 



(r(t))=0, (T(t)T(s))=S(t-s). 



(5) 



According to the dissipation-fluctuation theorem, the ther- 
mal noise intensity D is related to the friction coefficient 
and temperature T of the system, i.e., 



D = jk B T, 



(6) 



where fee stands for the Boltzmann constant. 

The external time-dependent force g(t) can be of any 
genre, both deterministic or stochastic (Luczka, 1999; Kula 
et al., 1998). This force makes the system to be in a non- 
equilibrium state and is a source of energy for movement of 
motors. For biological motors it comes from chemical reac- 
tions. Here, as an example and for the didactic purposes, 
we choose a time periodic force, namely 



g(t) = Acos(Clt), 



(7) 



where A is the amplitude and stands for the frequency of 
the external stimulus. Equivalently we can define the period 
T = 2-k/Q, of the external force. This kind of force can 
be realized for artificial motors while for biological motors, 
random non-equilibrium force is more adequate. 

2.1. Time scales 

In physics only relations between scales of length, time 
and energy are relevant, not their absolute values. There- 
fore we shell now translate the above equation of motion 
into its dimensionless form. First of all we determine char- 
acteristic quantities - time and length. The characteristic 
length is the period L of the potential V(x) and accordingly 
the coordinate of the molecular motor can be scaled as 



V = 



(8) 



Time can be scaled in several ways. One of the possibili- 
ties is the relaxation time tl of velocity or the correlation 
time of the velocity of the massive Brownian particle mov- 
ing freely and driven only by the thermal noise. It can be 
extracted from a special case of Eq. (1) when the right hand 
side is zero, i.e., 



mv + 77; = 0. 

From the above equation it follows 
v(t)=v(0) cxp(-i/T L ), 



(9) 



(10) 



where the characteristic time 77, = rn/7 is sometimes called 
the Langevin time. Another characteristic time comes from 
the overdamped motion of the particle in the potential 
V(x), when Eq. (1) reduces to the form 



dx 
7 ~dt = 



dVjx) 
dx 



(11) 



Then, by inserting into above equation the characteristic 
quantities, we get the definition of time tq, 



L AV 

7— = 

T 



L 



TO 



111 
AV" 



(12) 



During the time interval To overdampcd particle proceeds 
the distance L under the influence of the constant force 
AV/L. Third characteristic time follows from the friction 
less equation of motion, i.e., when Eq. (1) takes the form 



m- 



dt 2 



dV(x) 
dx 



(13) 



From the above equation the characteristic time r TO is given 
by the relations 



L AV 

tL L 



T 2 = 



mL 2 
AV 



(14) 



One can distinguish also other characteristic times like time 
period of the external driving T or the well known Einstein 
diffusion time 



TE = 



2D f 



D K = 



k B T 

7 



(15) 



Now we can rescale the equation of motion for the massive 
Brownian particle in several ways. Doing this we shell take 
as the relevant time scales these times which differs for dif- 
ferent systems, like r m or tq. Let us note that the Langevin 
time and the Einstein diffusion time do not depend on the 
system itself, i.e., on the potential and the external driving 
forces. 

2.2. Rescaled equations of motion 

Let us first propose to as the relevant time scale, i.e. the 
rescaled dimensionless time is s — t/r n . The dimcnsionless 
form of the Langevin equation then reads 



ey(s) + y(s) = F(y) + G(s) + ^/2D ^s), 



(16) 



where dot denotes derivative with respect to the rescaled 
times, F(y) = —dW(y)/dy — —W'(y) denotes the rescaled 
potential force and G(s) = (L/AV)g(t) = acos(ws) stands 
for the rescaled external driving force with the rescaled 
amplitude a = (L/AV)A and frequency uj = S1t . The 
rescaled spatially periodic potential W(y) = V(x)/AV = 
V(Ly)/AV — W(y + 1) has the unit period and the unit 
barrier height. The dimensionless mass 



m 
7 T o 



T_L_ 
T 



(17) 



is a ratio of the two characteristic times. The rescaled ther- 
mal noise £(s) = (L/AV)T(t) = (L/AV)T(t s) has exactly 
the same statistical properties as T(t). The dimensionless 
noise intensity 



D a = 



AV 



(18) 



is a ratio of thermal energy to activation energy the particle 
needs to traverse the non-rescaled potential barrier. 

On the other hand, if we choose the dimcnsionless time 
u = t/i~ m and r m as a time scale then we end up with 
another version of the rescaled Langevin equation, namely 



y(u) + jy(u) = F(y) + G(u) + ^2jD ^(u). 



(19) 



The dimensionless friction coefficient is a ratio of two char- 
acteristic times, different then previous two, namely, 



- T m T m 

m t l 



(20) 



The rescaled driving force G(u) = acos(u;u) with the 
rescaled amplitude a = (L/AV)A and frequency lj = f2r m . 

Two different scaling are useful in two limiting regimes: 
Eq. (16) in the overdamped case (when e « 1 is a small 
parameter) while Eq. (19) in the underdamped case (when 
-y << 1 is a small parameter). 

2.3. Numerical values for kinesin 

Let us evaluate characteristic times for one of the best 
known biological motor, namely, kinesin which moves along 
a microtubule. As already mentioned in the Introduction, 
microtubules are spatially periodic structures of period L w 
8nm. The mass of the kinesin head domain is of order m = 
lOOfcDa = 1.66 • \Q- 22 kg and its radius is R = 3nm. The 
friction coefficient 7 = 6- lCP n fcg/s is calculated from 
the Stokes formula with the use of the viscosity of water 
(77 = 1CT 3 kg /ms). In a typical Brownian domain the ac- 
tivation energy is 5 time higher than the thermal energy, 
AV = 5ksT and the temperature inside cell is about 310-ftT 
(37°C) . W can now estimate the typical characteristic times 
for the kinesin moving inside human cell; thus 

t l = 2.77- l(T 12 s, t = 1.8 • l(T 7 s, 
T m = 7 ■ l(T 10 s, t e = 4.57 • l(T 7 s. 

The dimensionless mass (17) in the first scaling and the 
dimensionless friction (20) in the second scaling have the 
following values 

s = 1.54 • 10" 5 < 1, 7 = 2.5-10 2 . 

One can note that the value of the parameter e is very 
small. Therefore Eq. (16) seems to be more appropriate 
than Eq. (19) because (16) contains the small parameter 
e. This allows, with a very good approximation, to put 
formally s — in the dimensionless equation of motion (16) 
yielding 



y(s) = F(y) + G(s) + y/2Dtf{s). 



(21) 



Analysis of this equation is much easier than Eq. (16) . From 
(21) it follows that the case of overdamped dynamics takes 
place for biological motors. On the other hand, Eq. (19) 
contains the parameter 7 and one cannot find justifiable 
arguments to neglect any term in it. So, this equation is 



Q 



much harder to analyze in a complete manner. However, 
if one wants to investigate some particular effects like the 
influence of the inertia on the transport, Eq. (19) is more 
practicable. For other motors, especially non-biological, in- 
ertial effect can be crucial. An example is an atomic Brow- 
nian motor moving in optical lattices (Brown and Renzoni, 
2008) and Eq. (1) in (Hagman ct al., 2008). It is worth 
to stress that if the non-equilibrium driving is of the form 
(7), the dynamics can be chaotic in some regimes leading 
to anomalous transport behavior like negative mobility or 
negative conductivity (Machura et al., 2007; Kostur et al., 
2008) . Below, we present characteristics of transport in full 
regime, from underdampcd to overdamped one. From now 
on we will use the later scaling (19) while defining perfor- 
mance characteristics of the motion of the Brownian mo- 
tor and shell omit all the hats in Eq. (19), for the sake of 
simplicity. 



3. Performance characteristics of molecular motor 

When one study the motion of molecular motors, the 
most important transport measure is an average stationary 
velocity (v) of the motor (Machura et al., 2004). Averaging 
should be performed over all realization of thermal noise, 
initial conditions and over a period of the external driving. 
Average velocity describes how much time a typical parti- 
cle needs to overcome a given distance in the asymptotic 
(long-time) state. This average velocity, however, is not the 
only transport attribute. Other characteristics are also im- 
portant. We will analyze two following transport aspects: 
quality and energetic efficiency. The quality of transport 
can be characterized by the effective diffusion coefficient 
D e jf, i.e., by the fluctuations in the position space (Lind- 
ner et al., 2001) 



Deff = 



2t 



(22) 



where the brackets (. . .) denote averaging over all realiza- 
tions of thermal noise and initial conditions. 

The third quantifier is the efficiency of noise rectification 
(Machura et al., 2004; Suzuki & Munakata, 2003). The mo- 
tor moves in viscous media. Therefore the minimal energy 
input required to move a particle in presence of friction 7 
over a given distance depends on the velocity. In this case, 
the rectification efficiency is given by the formula (Machura 
et al., 2005; Suzuki & Munakata, 2003; Linke et al., 2005; 
Kostur et al., 2006) 



V 



(23) 



The average input power P in corresponding to system de- 
scribed by Eq. (19) is given by the formula 



Pin - j[(v 2 ) ~ D ], 



(24) 



where D is defined in Eq. (18). This expression follows 
from the energy balance of the underlying equation of mo- 
tion (19) (Machura et al., 2004). 



4. Numerical experiment and discussion 

To examplify the above ideas we analyze Eq. (19) in the 
long-time asymptotic, time -periodic regime after effects of 
the initial conditions and transient processes have died out. 
Then, the statistical quantifiers of interest can be deter- 
mined in terms of the statistical average over the different 
realizations of the process (19) and over the driving period 
T. Since there are no analytical methods to handle (19), we 
made use of numerical approach and have carried out exten- 
sive and precise numerical simulations, applying Stochastic 
Runge-Kutta algorithm of order 2. For the initial condi- 
tions we took points from uniformly distributed circle in 
the phase space (x,v) with the radius r = \J x 1 + v 2 = 1 
and the origin in the point (x min , 0), where x m i n denotes 
minimum of the potential W(y). For the illustration of the 
above idea we fixed the following set of parameters: u> = 
4.9, D = 0.001,0.05 and the potential profile W{y) = 
AW[sin(27ry) + 0.25 sin(47ry)], where AW = 0.454 reduces 
the maximal barrier height to unity. The forces correspond- 
ing to this potential ranges from the minimal value -2.14 to 
the maximal value 4.28. If the amplitude a of the driving is 
higher than 4.28, the motor is able to overcome the poten- 
tial barrier in any direction just with the use of the external 
driving force and without of thermal noise. We analyze the 
system (19) in the 2-dimensional parameter space {7, a}. 
We change the rescaled friction coefficient 7 from 7 = 0.001 
(underdamped dynamics) to 7 = 100 (overdamped dynam- 
ics). The external time periodic force strength a varies be- 
tween and 5, where the later is just above the maximal 
value of the potential force. 

In Fig. 1 we show three characteristics in the low tem- 
perature regime, D = 0.001, cf. Eq. (18). In panel (a), we 
depict the stationary average velocity of the motor. White 
area corresponds to the case when the absolute value of the 
velocity is small, |(i>)| < 10~ 2 . Blue (or gray) areas sur- 
rounded by contours stands for the negative mean velocity 
and red (gray) color indicates regions of positive velocity. 
It is seen that for a chosen set of parameters the Brownian 
motor has a noticeable velocity for the force amplitude a > 
2 and for not strong damping 7 < 2. Strictly speaking, the 
velocity is not exactly zero assuming very small values of 
order 10~ 2 -j- 10~ 7 (white area on Fig. 1 (a)). By inspecting 
the several colored (gray) areas, one can observe that the 
velocity as a function of 7 or the amplitude a exhibits the 
multiple velocity reversals (Kostur and Luczka, 2001). Be- 
cause the friction coefficient 7 depends on the linear size of 
the motor, it means that motors of different sizes can move 
in opposite directions. In panel (b) we show the rectification 
efficiency (23). There, white area corresponds to efficiency 
smaller that in < 10 -3 . The correlation of dark regions in 
panels (a) and (b) is evident. There are two main islands of 
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Fig. 1. (color online) Performance characteristics of the molecular 
motor are presented versus the external force amplitude a and friction 
coefficient 7. In panel (a) we show the average velocity \{v)\ < 10 — 2 
(white), negative (blue, gray surrounded by contour), positive (red, 
gray). In panel (b) we present the efficiency. In panel (c) we depict 
the effective diffusion. Other parameters are: Do = 0.001, ui = 4.9. 

efficient energy conversion - one for {7, a} £ {0.4 2, 3^-5} 
and second for {7, a} £ {0.001 -j- 0.2, 2 5}. In the region 
of parameters {7 « 0.002, a w 2.9}, the motor efficiency 
is very high and is almost 90% (remember that the sys- 
tem is far from an equilibrium state) . In panel (c) we show 
the effective diffusion coefficient of the motor. Large value 
of this quantifier means low quality of motion - Brownian 
particles move in a very irregular manner. This region of 
high efficiency and low quality of motion (high diffusion), 
makes our stochastic model similar to deterministic (Ma- 
teos, 2000) or Hamiltonian (Schanz et al., 2001) systems. 

Upon the inspection of Fig. 1 one can notice that for the 
parameter island with large negative average velocity and 
high efficiency (large area around 7 = 0.02 and a = 2.8), 
the motor moves quite irregularly because -D e // is large. 
For the island with large positive velocity and medium effi- 
ciency (smaller area around 7 = 1 and a = 4), the effective 
diffusion is small reflecting good quality of transport and 
regular motion. 

Let us consider the case of higher temperature, D Q = 
0.05. Transitions of a particle over potential barriers are 
easier and details of the shape of the potential is not so 
important now. As one can see from panel (a) of Fig. 2, ab- 
solute values of the average velocity are one order smaller 



Fig. 2. (color online) The same as in Fig. 1 but for higher temperature, 
Dq = 0.05. Panel (a): the average velocity \(v)\ < 10~ 3 (white), 
negative (blue, gray surrounded by contour), positive (red, gray). 
Panel (b): the efficiency. Panel (c): the effective diffusion. 



than in the previous temperature regime. Higher temper- 
ature makes the motor motion slower in both negative - 
blue (gray) contoured areas, or positive - red (gray) col- 
ored regions. Again as in the previous case, in the the 2- 
dimensional parameter space {7, a}, we can identify the ve- 
locity reversal phenomenon upon the change of the control 
parameters a or 7. In the negative direction the maximal ve- 
locity (v) — —0.06 and in the positive direction (v) = 0.14. 
In the optimal regime, the efficiency of the motor transport 
is of order 1%. The only characteristic which seems to be 
better than in the previous scenario is the effective diffu- 
sion presented in panel (c). The motor seems to move in 
a much more regular manner than for lower temperature 
yielding maximum value of -D e // = 250 for underdamped 
system. It is a region where the averaged velocity is almost 
zero and the efficiency is very low. 

Comparison of Figs 1 and 2 leads to the conclusion that 
the influence of higher temperature is rather destructive. 
The only better quantifier is the effective diffusion coef- 
ficient which is smaller. It means that the islands in the 
parameter space {7, a} for higher temperature with rela- 
tively high (negative or positive) velocity accompanied by 
the highest possible efficiency represent regimes where the 
motor moves rather in a regular way. 



r. 



Conclusions 

The most demanded properties of any transporting ma- 
chinery are: efficiency and quality. The latter in our case 
is characterized by the effective diffusion coefficient. Let 
us note that even in any equilibrium system the particle 
can be transported over long distances due to thermal dif- 
fusion. It is, however, very unreliable when the distance 
becomes large because the diffusion cannot distinguish di- 
rection and most of traveling particles would not arrive at 
prescribed destination in a reasonable time. On the other 
hand, from Figs 1 and 2 we can see that the large diffusion 
is sometimes in regimes where the energetic efficiency is 
high. In real situations, Nature chooses between above sce- 
narios. Under some circumstances, if e.g. the distances are 
small, the diffusion can be exploited to transport the cargo 
and no energy is wasted. When the distance is larger, then 
the non-equilibrium transport is applied and reliability is 
achieved sacrificing the energy input. 
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